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Lattice QCD allows computations of moments of structure functions 
from first principles. An overview of the present status of the calculations 
is given. Recent results and future perspectives are discussed. 
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1. Introduction 

The moments of structure functions can be calculated from first prin- 
ciples, without making any model assumptions, using the non-perturbative 
techniques of lattice QCD. Numerous results have been obtained from lattice 
computations in the last years. They include the calculation of the lowest 
moments of the unpolarized structure functions, of the spin-dependent gi 
and g2 structure functions and of the h\ transversity structure function. 
Some higher-twist matrix elements have been studied as well. 

The novelties of the last couple of years are the first computations of 
structure functions done in full QCD (unquenched), and new proposals for 
the extrapolations to the chiral limit using chiral perturbation theory. A lot 
of work however still needs to be done in order to control other systematic 
uncertainties like the continuum limit or the non-perturbative renormaliza- 
tion. 



2. Structure functions on the lattice 

It is not possible to compute the structure functions directly on the 
lattice, as they describe the physics close to the light cone, and Monte 
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Carlo simulations are instead done in Euclidean space. One can however 
calculate their moments via OPEs which have the general form 

C x n Fi(x,Q 2 ) ~ C\ n (Q 2 /^ 2 ) ■ (h\OM\h). (1) 
Jo 

The Wilson coefficients Cj jfl contain the short-distance physics, calculable 
in continuum perturbation theory. The long-distance physics is contained 
in the hadronic matrix elements, which can be computed doing Monte Carlo 
simulations and then performing a proper lattice renormalization. 

The expansions of Eq. (1) are dominated by operators of twist two. 
The moments of the unpolarized quark distributions, of the spin-dependent 
structure function g\ , of the component of <?2 and of the h\ transversity 
structure function are respectively measured by the following operators: 

(x n ) ~ <fc|^7{A-" £) M^I /l >> ( 2a ) 

((Ax)") ~ (h\^ l5l{lx D fll ---D^ n} ^\h), (2b) 

(x n ) d2 ~ (h^-f^D^D^yD^y^lh), (2c) 

((Sx) n ) ~ (h^jsa^D^-.-D^^lh). (2d) 

The moments of the unpolarized gluon distribution, 

(x n ) g ~ ( h | ]T Tr • • • D^Fp^) \ h), (3) 

p 

which in the flavor singlet case mix with the moments in Eq. (2a), have 
been quite hard to compute on the lattice up to now, and for this reason 
only non-singlet quark distributions are generally considered. 

Structure functions calculations on the lattice have been done with Wil- 
son fermions, which however break chiral symmetry. This causes an additive 
mass renormalization even for vanishing bare masses, and a heavy lattice 
pion (m n ~ 500 MeV). Therefore one needs to perform extrapolations to 
the chiral limit, and this systematic error has to be controlled. 1 

Extrapolations to the continuum limit a — ► need also to be done. To 
get a faster convergence to the continuum and decrease the systematic error 
arising from the finiteness of the lattice spacing, 0(a) improvement has been 
implemented in most cases: the contributions of 0(a) are removed so that 
one has 
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1 Recent formulations of chiral fermions on the lattice, called Ginsparg- Wilson (in 
particular: overlap, domain- wall and fixed-point fermions), possess exact chiral sym- 
metry also at finite lattice spacing and thus do not suffer from this problem. 
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This is achieved in on-shell matrix elements by adding a counterterm to the 
Wilson action, AS{ = c sw ■ '\g a A r/4,a Y. x ^v ^(x)a flu F^' VCT (x)'ip(x). The 
improvement coefficient c sw , which has to be exactly tuned so as to cancel 
all 0(a) contributions coming from the action, is now well known. The var- 
ious operators have to be improved as well, and this is achieved by adding 
as counterterms bases of higher-dimensional operators with the same sym- 
metries as the original unimproved ones: O imp = (1 + bo am) O + a J2i c i Oi 
(with dim(Oj) = dim(O) + 1), and determining the coefficients bo and Cj. 

After the matrix elements have been simulated with Monte Carlo algo- 
rithms, they need to be renormalized from the lattice to the continuum: 

(or nt ) = E (** - (^'f - R T)) ■ (of')- (s) 

Since Lorentz and (for Wilson fermions) chiral symmetry are broken on the 
lattice, mixing of operators under renormalization is more complicated than 
in the continuum. Much effort has been spent on computing perturbative 
renormalization factors, which are now known for the lowest three moments 
of all structure functions in the case of Wilson (often in the 0(a) improved 
theory) and of overlap fermions. For renormalization calculations see Ref. [1] 
and references therein. The relevant perturbative renormalization factors 
have been calculated also for some classes of higher-twist operators [2]. 

Non-perturbative renormalization has instead been quite hard to do in 
the standard Wilson case, where it is difficult to find a plateau for extracting 
the signal. This is not the case in the Schrodinger Functional approach, 
where one has a better control over many systematic errors, although the 
computation of the running of renormalization factors is quite tedious [3] . 

For reasons of computing power, structure functions have been studied 
mostly in the quenched approximation, and only recently full QCD com- 
putations have become feasible. On the lattice the (Grassmann) fermion 
variables are integrated out, and in full QCD one uses the partition func- 
tion 

Z = J VU det(P[U] + m q ) e~ s ^ (U^ = e iagoA ^). (6) 

Quenching instead amounts to putting det(]ft[U] + m q ) = 1 in the simu- 
lations. In physical terms, this means that there are no sea quarks in the 
calculations: the internal quark loops are neglected. Although it looks quite 
drastic, in many cases this does not turn out to be a bad approximation. 

3. Results 

Various collaborations have been studying structure functions on the 
lattice in the last years: the QCDSF Collaboration, the LHPC Collabo- 



4 



CAPITANI PRINTED ON FEBRUARY 1, 2008 



ration, and a collaboration which uses the Schrodinger Functional scheme. 
See Refs. [4, 5, 3] respectively for recent works and further references. 

The major advancement in the last couple of years has been the compu- 
tation of the structure functions in full QCD. The first studies with dynam- 
ical fermions have shown that for various structure functions of the u and 
d quarks there are no statistically significant differences between quenched 
and full QCD results. This was discovered by the LHPC Collaboration and 
then confirmed by the QCDSF Collaboration. We then have two calcula- 
tions in full QCD with the same message: the results for quenched and 
unquenched moments are nearly equal, and the observed discrepancies of 
lattice results with experiment are still there. Quenching had been previ- 
ously conjectured as the cause of these discrepancies, but evidently this is 
not the case, at least for the values of quark masses attainable at present, 
which are unfortunately quite large. 

Usually only the lowest three or so moments of the various structure 
functions can be calculated, due to increasing computational difficulties and 
mixing problems. The largest discrepancies concern unpolarized quark dis- 
tributions. The spin-dependent g\ and g% structure functions have been 
extensively studied as well, and also the axial charge qa = Au — Ad, whose 
result is not so far from experiment. The gi structure function with Wilson 
fermions presents a mixing problem due to chirality breaking: 2 there is 
a mixing with lower-dimensional operators, which gives rise to power di- 
vergences ~ l/a n in the continuum limit. The hi transversity structure 
function, whose lowest moment is the tensor charge, has also been studied 
on the lattice. A recent lattice result for this quantity, which has never 
been measured, is Su — 5d = 1.21(4) [6]. An experimental measure of the 
hi transversity structure functions would therefore be quite interesting. 

The results coming from lattice QCD are getting more and more precise, 
so it would also be useful to have more precise measurements of parton 
distributions and above all a careful analysis of their errors. Some first 
studies in this direction have been presented at this conference [7]. 

Nevertheless, on current lattices there are still some practical limitations. 
One of them is given by the difficulty to compute disconnected diagrams 
(i.e. connected, but only by gluon lines). Only exploratory studies have 
been made so far, and these diagrams have not yet been included in lattice 
computations. They are however flavor independent, and therefore they do 
not contribute to the differences between u and d structure functions (for 
iTiu = iTid, which is the case). This means that quantities like qa = Au — Ad 
and {x n ) u -d do not receive contributions from disconnected diagrams. 

Another limitation is given by the extrapolations to the chiral limit. 



2 With Ginsparg- Wilson fermions this mixing is thus forbidden. 
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Chiral and continuum extrapolations are performed using the fit formula 

A + Bml + Ca 2 , (7) 

and doing just the naive chiral extrapolations in m\ could also turn out to be 
an explanation for the discrepancies. The extrapolations to the continuum 
limit are quadratic in the lattice spacing when improved fermions are used. 

There have also been quenched calculations which use the Schrodinger 
Functional. This is a finite volume scheme (with Dirichlet boundary con- 
ditions on the time direction) in which it is possible to do simulations at 
very small quark masses, and the pion is much lighter than in the stan- 
dard Wilson case. The renormalization scale is identified with the inverse 
size of the lattice, fi = 1/L, which allows one to apply finite-size scal- 
ing techniques; there is no need to go to the infinite volume limit. In the 
Schrodinger Functional calculations the renormalization of operators can be 
done non-perturbatively, at scales extending over some orders of magnitude. 
The results for the lowest moment of the unpolarized structure function of 
the pion, in the MS scheme at the renormalization scale \i = 2.4 GeV, is 
(x) = 0.30 ± 0.03, while the experimental number is (x) = 0.23 ± 0.02. 



4. Chiral extrapolations 



It has recently been suggested that extrapolations of the lattice data us- 
ing chiral perturbation theory could solve the discrepancy with experiment. 

It seems that the pion cloud of the nucleon is not adequately described 
by current lattices. The IR behavior of pion loops generates chiral logs 
~ m\ lnm 2 , and using chiral perturbation theory [8] one gets for the unpo- 



larized case (x r 



A. 



Introducing 



1 - (3g\ + l)/(4vr/ u ) 2 • m\ lnm 2 

in this non-analytic term a phenomenological cutoff A (the size of the source 
generating the pion cloud) one then arrives at the fit formula [9] 



u—d 



(ggi + 1) 

(4tt/ w )2 



ml In 



ml 



m 2 + A 2 



+ B n ml + C n a 2 . (8) 



It looks like phenomenological extrapolations which use chiral pertur- 
bation theory could solve the discrepancy with experiment. However, they 
need to be better investigated, and although they have the potential to 
explain the discrepancy, it is still too soon to say something certain. 3 

So far we are able to reproduce the phenomenological results only at 
posteriori. The corresponding values of A lie between 300 MeV and 700 



3 A recent calculation [10] seems also to show that in the polarized case the next order 
corrections give a fit formula which deviates very little from the naive extrapolations. 
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MeV, so A does not seem at present to have some predictive power. More- 
over, the currently available lattice data do not even discriminate between 
naive chiral fits, Eq. (7), and chiral perturbation theory fits, Eq. (8). The 
problem is that in present simulations the pions are not sufficiently light. 
A smaller pion mass is needed (m n < 250 MeV) in order that the param- 
eters of the chiral expansions can be well determined on the lattice. For 
such determinations a computing power of about 8 Teraflops for one year is 
required [9]. The next generation of computers, which will come in a couple 
of years, should then be able to perform these calculations. 

For the pion cloud of the proton to be adequately included in the lattice 
box and properly measured in Monte Carlo simulations, the pion correlation 
length should also be smaller than the lattice size, otherwise it will not be 
fully contained in the lattice. The length of a spatial dimension for currently 
available lattices is about 2 fermi, so it would be very useful to do simulations 
also with larger physical volumes and minimize finite volume effects. 4 

The pion cloud will be the focus of many future lattice investigations. 



5. Higher twist 

There are a few lattice results regarding higher-twist corrections for the 
pion and the nucleon [2]. They show that, for the twist- four contributions 
coming from 4-fermion operators like J2a ^7^75 t Ar( P ' V^75 t A ip, the 1/Q 2 
power corrections are, at least in these particular cases, quite small. 

The calculations give, for a particular contribution (I = 2) to the twist-4 
matrix element of the first moment of the unpolarized pion distribution, 

1 dxF 2 (x,Q 2 ) = 1.67(64) ■ a s {Q 2 ) ■ /f, (9) 
io 1=2 Q 

and for the proton, for a particular class of operators (flavor 27, 1=1), 



C dxF 2 (x,Q 2 ) = -0.0006(5) • « 5 (Q 2 ) • (10) 

J0 27,7=1 y 

These special flavor and isospin combinations are chosen to avoid mixing 
with lower-dimensional operators. The above twist-4 numbers are much 
smaller than the twist-2 lattice matrix element (at Q 2 = 4 GeV 2 ) of the 
pion, Jg dx F2(x,Q 2 ) = 0.152(7), and of the proton, which is about 0.14. 

We point out that these studies, contrary to leading-twist operators, 
are not complete and systematic. Mixing problems have so far limited the 
calculations only to a restricted class of operators. 

4 We should also mention that the lattice results for elastic form factors seem to have 
discrepancy problems similar to the structure functions, particularly in the case of 
the magnetic form factor. It is likely that in this case too the pion cloud is not being 
properly simulated and finite-size effects are present. 
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6. Perspectives 

The lattice calculations of structure functions are getting more and more 
refined. However, the physics of the pion cloud must be properly included, 
and simulations on lattices of large physical volumes as well as proper ex- 
trapolations to the chiral limit will need to be done. 

Ginsparg- Wilson fermions possess exact chiral symmetry and could help 
to study the chiral limit. Calculations with the Schrodinger Functional will 
also be useful for a better understanding of the approach to the chiral limit. 

Contributions of higher-twist operators to moments of parton distribu- 
tions seem to be quite small, but they have not been systematically studied. 
They will continue to be a challenge for lattice QCD for a few more years. 

Disconnected diagrams will have to be included in the simulations. 

An almost virgin territory for the lattice is the study of the momentum 
and spin distributions of the gluon, which so far have been plagued in Monte 
Carlo simulations by large statistical errors. It would be quite interesting 
to know them. The lattice distributions of sea quarks are also unknown. 

REFERENCES 

[1] S. Capitani, G. Rossi, Nucl. Phys. B433, 351 (1995); 

G. Beccarini, M. Bianchi, S. Capitani, G. Rossi, Nucl. Phys. B456, 271 (1995); 

M. Gockeler et al., Nucl. Phys. B472, 309 (1996); 

S. Capitani, Nucl. Phys. B592, 183 (2001); 

S. Capitani, Nucl. Phys. B597, 313 (2001). 
[2] S. Capitani et al. (QCDSF), Nucl. Phys. B570, 393 (2000); 

S. Capitani et al. (QCDSF), Nucl. Phys. B (Proc. Suppl.) 94, 299 (2001). 
[3] M. Guagnelli, K. Jansen, R. Petronzio, Phys. Lett. B493, 77 (2000); 

K. Jansen, hcp-lat/0010038 (ICHEP 2000, Osaka). 
[4] M. Gockeler et al. (QCDSF), hep-ph/0108105; 

S. Capitani et al. (QCDSF), Nucl. Phys. B (Proc. Suppl.) 106, 299 (2002); 

see also: C. Best et al. (QCDSF), Phys. Rev. D56, 2743 (1997). 
[5] D. Dolgov et al. (LHPC), Nucl. Phys. B (Proc. Suppl.) 94, 303 (2001); 

D. Dolgov et al. (LHPC), hep-lat/0201021, to appear in Phys. Rev. D. 
[6] S. Capitani et al. (QCDSF), Nucl. Phys. B (Proc. Suppl.) 79, 548 (1999). 
[7] In these proceedings: J. Blumlein; H. Bottcher; J. Soffer; R.S. Thornc; 

see also: R.S. Thorne, hep-ph/0205235; R.S. Thorne et al., hep-ph/0205233; 

J. Blumlein, H. Bottcher, hep-ph/0203155, to appear in Nucl. Phys. B. 
[8] A.W. Thomas, W. Melnitchouk, F.M. Steffens, Phys. Rev. Lett. 85, 2892 

(2000); D. Arndt, M.J. Savage, Nucl. Phys. A697, 429 (2002); J.-W. Chen, 

X. Ji, Phys. Lett. B523, 107 (2001); Phys. Rev. Lett. 88, 052003 (2002). 
[9] W. Dctmold et al., Phys. Rev. Lett. 87, 172001 (2001). 
[10] W. Dctmold, W. Melnitchouk, A.W. Thomas, hcp-lat/0206001. 



